Your thesaurus codes are: 
(06.05.1; 06.07.1; 06.16.1) 



ASTROPHYSICS 
25.5.1996 



Updating Standard Solar Models 

F. Ciacio\-^ S. DegrinnocentP,^ and B. Ricci^ 

^ Dipartimento di Fisica deU'Universita di Ferrara, via Paradise 12, 1-44100 Ferrara, Italy 

^ Max-Planck Institut for Astrophysics, K. Schwarzschild Str. 1, D-85740 Garching bei Miinchen, Germany 

^ Istituto Nazionale di Fisica Nucleare, Sezione di Ferrara, via Paradiso 12, 1-44100 Ferrara, Italy 

Received ; accepted 



Abstract. We present an updated version of our standard 
solar model (SSM) where helium and heavy elements diffu- 
sion is included and the improved OPAL equation of state 
(Rogers 1994, Rogers Swenson & Iglesias 1996) is used. In 
such a way the EOS is consistent with the adopted opac- 
ity tables, from the same Livermore group, an occurrence 

0^ which should further enhance the reliability of the model. 

On The results for the physical characteristics and the neu- 
trino production of our SSM are discussed and compared 
^^with previous works on the matter. 

Key words: The Sun: evolution - The Sun: general - The 
Sun: particle emission 
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In the last decades evolutionary computations of Stan- 
T— I dard Solar Models (SSM) played a fundamental role in un- 
*0 derstanding the inner solar structure and, in particular, 



In addition, updated values of the relevant nuclear cross 
sections are used and more refined values of the solar con- 
stant and age (BP95) are adopted. 

Before entering into the argument, let us recall that in 
recent years helioseismology has added important pieces 
of information on the solar structure, producing severe 
tests for standard solar model calculations. According to 
Christensen-Dalsgaard et al. (1993) one can accurately de- 
termine the depth of the solar convective zone and the 
speed of sound at its bottom: 



Rt/Re = 0.710^0.716 

Cb = 0.221^ 0.225 (Mm/s) 



(1) 



in approaching the well known problem of solar neutri- 



nos. Bahcall & Pinsonneault (1995; hereafter BP95) have 
already shown that to reach the agreement with obser- 



Richard et al. (1996; hereafter RCVD96) recently con- 
firmed the value of R;,, finding R;,/R0=O.7137. In addition 
several determinations of the helium photospheric abun- 
dances have been derived from inversion of helioseismo- 
logical data, with results in the range (see Castehani et 
al. 1996 and refs. therein): 



Oh, 

o 



vational evidence given by helioseismology one needs, in ^photo — 0.233 . 0.268. 



(2) 



'addition to the best available input physics, an accurate 
treatment of element diffusion all along the solar struc- 
ture. The recent availability of an improved equation of 
state (EOS), as given by the OPAL group (Rogers 1994, 
Rogers Swenson & Iglesias 1996) obviously suggests to in- 
vestigate the influence of such an improvement on SSM. In 
this paper we discuss this scenario, presenting an updated 
SSM resulting from a recent version of FRANEC (Frascati 
Raphson Newton Evolutionary Code), where helium and 
heavy elements diffusion is included and the OPAL equa- 
tion of state is used. Note that in such a way the EOS is 
consistent with the adopted opacity tables, from the same 
Livermore group (Rogers & Iglesias 1995), an occurrence 
which should further enhance the reliability of the model. 
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The much smaller errors often quoted, reflect the observa- 
tional frequency errors only. The results actually depend 
on the method of inversion and on the starting physical 
inputs (e.g. the EOS), see RCVD96. 

Note that, since in building a SSM one is dealing with 
three free parameters (mixing length, original He content 
and original Z/X), with these three additional constraints 
(R;,, c;, and ^ photo) no free parameter is left for SSM 
builders. 

After discussing the effect of the physical inputs, we 
compare our SSM with other recent solar model calcu- 
lations, all including diffusion of helium and heavy ele- 
ments, finding an excellent agreement. We present neu- 
trino fiuxes and the expected signals in ongoing ex- 
periments. A detailed presentation of our new model 



tion...) is available on World Wide Web at the address 
http://dns.unife.it/Fiorentini/index.htm and a more ex- 
tensive discussion will be published elsewhere. 

As for tho computations, FRANEC has been described 
in previous papers (see e.g. CliiefR & Straniero 1989, 
Castellani et al. 1992). Recent determinations of the the 
solar luminosity {Lq = 3.844 • 10^^ erg/s) (BP95) and 
of the solar age (tQ = 4.57 • 10'' yr) (BP95) are used. 
The present ratio of the solar metallicity to solar hydro- 
gen abundance by mass corresponds to the most recent 
value by Grevesse & Noels (1993) : (Z/X)photo= 0.0245. 

Following the standard procedure (see e.g. Bahcall & 
Ulrich 1988), for each set of assumed physical inputs the 
initial Y, Z and the mixing length parameter a were var- 
ied until the radius, luminosity and Z/X at the solar age 
matched the observed values within a tenth of percent or 
better. We considered the following steps: 

a) As a starting model we used the Straniero (1988) 
equation of state, the version of the OPAL opacity tables 
available in 1993 (Rogers & Iglesias 1992, Iglesias Rogers 
& Wilson 1992) for the Grevesse & Noels (1993) solar 
metallicity ratio, combined with the molecular opacities 
by Alexander & Fergusson (1994); diffusion was ignored. 
This model is useful for a comparison with BP95 (without 
diffusion) which uses the same chemical composition. 

b) Next, we introduced the OPAL equation of state. 
With respect to other commonly used EOS, this one 
avoids an ad hoc treatment of the pressure ionization 
and it provides a systematic expansion in the Coulomb 
coupling parameter that includes various quantum effects 
generally not included in other computations (see Rogers 
1994, Rogers Swenson & Iglesias 1996 for more details). 

c) Next, OPAL 1993 tables were substituted with the 
latest OPAL opacity tables (Rogers & Iglesias 1995), again 
for Grevesse & Noels (1993) solar metallicity ratio. With 
respect to Rogers & Iglesias (1992) the new OPAL ta- 
bles include the effects on the opacity of seven additional 
elements and some minor physics changes; moreover the 
temperature grid has been made denser. 

d) Furthermore, we included the diffusion of helium 
and heavy elements. The diffusion coefficients have been 
calculated using the subroutine developed by Thoul (see 
Thoul, Bahcall and Loeb 1994). The variations of the 
abundance of H, He , C, N, O and Fe are fohowed ah along 
the solar structure; ah these elements are treated as fuhy 
ionized. According to Thoul, Bahcall and Loeb (1994) all 
other elements are assumed to diffuse at the same rate as 
the fully-ionized iron. To account for the effect of heavy el- 
ement diffusion on the opacity coefficients we interpolated 



with different total metallicity (Z=0.01,0.02,0.03,0.04). 

e) As a final point, we investigated the effect of updat- 
ing the nuclear cross sections for ^He-|-^He and ^He-|-^He 
reactions, following a recent reanalysis of all available data 
(see Castehani et al. 1996, and table 6). For 5(0)34, a poly- 
nomial fit gives (energy in Mev, S in Mev barn): 

SaiiE) = (4.8 - 2.9E + 0.9^^) x 10"^ (3) 

Alternatively, by using an exponential parametrization, 
as frequently adopted in the literature, one obtains an 
equally good fit to experimental data: 

SsiiE) = 5.1 ■ 10-^ea;p(-0.83£ -I- 0.25£^). (4) 

At the energies of interest, Eq « 20 Kev, the second ex- 
pression yields an S factor larger by about 5%, a value 
which is indicative of the uncertainty on S34 at these low 
energies. In the calculations, we used eq. 4. 

The resulting solar models are summarized in table 1, 
which deserves the following comments: 

(a — )■ b): the introduction of the new OPAL EOS 
reduces appreciably the initial helium abundance. The 
Straniero (1988) EOS understimates the Coulomb effects 
neglecting the contribution due to the electrons, which are 
considered as completely degenerate, whereas the OPAL 
EOS includes corrections for Coulomb forces which are 
correctly treated (see Rogers 1994 for more details). The 
models with Straniero EOS have a higher central pressure 
and a higher central temperature, and correspondingly a 
higher initial helium abundance. The effect of an understi- 
mated Coulomb correction was discussed in Turck-Chieze 
& Lopez (1993). Note that the hehum abundance is higher 
than the helioseismological determination by more than 
2a. Moreover, the transition between radiative and con- 
vective regions is not correctly predicted by the model, 
the convective region being definitely too shallow. 

(b c) : the updating of the radiative opacity coeffi- 
cients has minor effects. The surface helium abundance is 
now within 2a from the helioseismological determination, 
but the convective zone is again too shallow. 

(c — )• d): this step shows the effects of diffusion. He- 
lium and heavy elements sink relative to hydrogen in the 
radiative interior of the star because of the combined ef- 
fect of gravitational settling and of thermal diffusion. This 
increases the molecular weight in the core and thus the 
central temperature is raised. The surface abundances of 
hydrogen, helium and heavy elements are appreciably af- 
fected by diffusion: the initial value Yj„= 0.269 is reduced 
to the present photospheric value Yp;jo(o=0.238, within 
2a from the hehoseismological results. More important. 



speed are now in good agreement with helioseismologi- 
cal values. With respect to models without diffusion, in 
the external regions the present helium fraction is re- 
duced while the metal fraction stays at the observed pho- 
tospheric value. Thus the opacity increases and convection 
starts deeper in the Sun. As an example of the diffusion 
process we show in fig.l the time dependence of the He 
profile and in fig. 2 the present H profile, calculated with 
and without diffusion. 

(d e): The modifications of our Solar Model aris- 
ing from the new values of the nuclear cross sections are 
negligible with respect to the other improvements just pre- 
sented. 



Fig. 1. Time dependence of the He profile of our model e). 
The He profile at t=4.57 Gyr is that of our "best" standard 
solar model 



Fig. 2. H profile for a solar model without diffusion (our model 
c) and with diffusion (our model e). 

Our "best" Standard Solar Model, model (e), appears 
in good agreement with recent calculations (table 3) by 
several authors, all including microscopic diffusion and us- 
ing (slightly) different physical and chemical inputs, sum- 
marized in table 4. Our ingredients are very close to those 
of BP95 and one notes a substantial similarity with that 
model. We only have a slightly lower central temperature, 
possibly an effect of the different EOS. Among the differ- 
ent models, diffusion looks less efficient in RVCD96, as a 
consequence of the inclusion of rotational induced mixing. 

The predicted neutrino fluxes and signals from our 
SSM are summarized in table 2. All in all, the results 
are quite stable with respect to the changements we have 
introduced as long as diffusion is neglected. On the other 
hand, due to the higher central temperature, model (d) 
has significantly higher and CNO neutrino fiuxes. It is 
essentially the increase of the boron fiux which enhances 
the predicted Chlorine signal. The slight changement in 
the nuclear cross section weakly affects neutrino fluxes and 
signals: the '^Be and ^B fluxes are reduced by about 5%, 



S34. Should we use the polinomial parametrization of eq. 
3 one would get a further 5% decrease. As a flnal remark, 
we want to stress that again our results are in excellent 
agreement with other recent calculations including diffu- 
sion (see table 5), but the prediction of Kovetz & Shaviv 
(1994) for the ^^N and ^^O fluxes which are an order of 
magnitude smaller than those of all other calculations for 
reasons which are unknown to us ^ . 

Regarding the comparison with BP95, the small dif- 
ference in the neutrino production reflects well the small 
difference in the central temperature, already remarked, 
and in the adopted values of S34 (see table 6) . 
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Table 1 : Comparison among solar models obtained with different versions of 
the FRANEC code. The labels (a) to (e) correspond to the models defined 
in the text. Our best Standard Solar Model is (e). The last column shows 
the helioseismological results as discussed in the text. Here: S88=Straniero 
(1988) 





(a) 


(b) 


(c) 


(d) 


(e) 
best 


Helioseism. 


i0[Gyr] 

L0[1O'' erg/cmVs] 
i?0[lO^° cm] 

a 


4.57 

3.846 

6.961 

0.0245 

2.023 


4.57 

3.843 

6.963 

0.0245 
1.774 


4.57 

3.844 

6.959 

0.0245 

1.786 


4.57 

3.843 

6.959 

0.0245 

1.904 


4.57 

3.844 

6.960 

0.0245 

1.901 




Y 

2 in 


0.699 

0.284 
0.0171 


0.718 

0.265 
0.0176 


0.722 
0.261 
0.0177 


0.711 

0.269 
0.0198 


0.711 
0.269 
0.0198 




^^photo 
^ photo 
'^photo 


0.699 
0.284 
0.0171 


0.718 
0.265 
0.0176 


0.722 
0.261 
0.0177 


0.743 
0.238 
0.0182 


0.744 
0.238 
0.0182 


0.233^ 0.268 


Rb/RQ 
T6[10'^ K] 
CblW'' cm s"^] 


0.738 

1.99 

2.11 


0.726 

2.10 

2.16 


0.728 

2.08 

2.16 


0.716 

2.17 

2.22 


0.716 

2.17 

2.22 


0.710^ 0.716 
2.21^ 2.25 


T,[10" K] 
p4100 gr cm" 3] 

Yc 


1.555 
1.524 
0.63 


1.545 
1.472 
0.61 


1.542 
1.470 
0.61 


1.569 
1.514 
0.63 


1.569 
1.518 
0.63 




OPACITY 
EOS 


OPAL 

S88 


OPAL 
OPAL 


OPAL 
OPAL 


OPAL 
OPAL 


OPAL 
OPAL 





Table 2: Neutrino fluxes and signals obtained with different versions of the 
FRANEC code. The labels (a) to (e) correspond to the models defined in 
the text. Our best prediction is (e). 











(a) 


(b) 


(c) 


(d) 


(e) 


















best 


$pp [10^ cm" 


-2 


S" 




60.17 


60.37 


60.66 


59.76 


59.92 


$pep[109 cm- 


-2 


S" 




0.14 


0.14 


0.14 


0.14 


0.14 


$Be [10^ cm- 


-2 


S" 




4.58 


4.22 


4.09 


4.71 


4.49 


$N [10^ cm" 


-2 


S" 




0.39 


0.36 


0.35 


0.52 


0.53 


$0 [10^ cm- 


-2 


S" 




0.33 


0.30 


0.29 


0.45 


0.45 


$B [10*^ cm" 


-2 


s" 




4.73 


4.18 


3.95 


5.37 


5.16 


SGa[SNU] 








126 


121 


120 


130 


128 


Sc;[SNU] 








6.9 


6.2 


5.9 


7.7 


7.4 



Table 3: Comparison among several SSMs, all including diffusion. The cor- 
respondence between acronyms and references is as follows: CGK89=Cox et 
al. (1989), P94=ProfRt (1994), KS94=Kovetz k Shaviv (1994), RVCD96= 
Richard et al. (1996). BP95 indicates here the "best model with diffusion" 
of Bahcall &: Pinsonneault 1995. FRANEC96 indicates our best model with 
diffusion, model (e). The star indicates values of cj calculated by us assum- 
ing fully ionized gas EOS. 





CGK89 


P94 


KS94 


RVCD96 


BP95 


FRANEC96 


<0[Gyr] 

Lq[W'' erg/cmVs] 
Rq[10^° cm] 

{Z/X)photo 


4.54 
3.828 
6.9599 
0.02464 


4.60 

3.846 

6.9599 

0.02694 


4.54 
3.8515 
6.960 
0.02763 


4.60 
3.851 
6.959 
0.0263 


4.57 
3.844 
6.9599 
0.02446 


4.57 
3.844 
6.960 
0.0245 


Yin 


0.691 
0.289 
0.02 


0.6984 
0.2803 
0.02127 


0.6797 
0.2991 
0.0211 


0.7012 
0.2793 
0.0195 


0.70247 
0.27753 
0.02 


0.711 
0.269 
0.0198 


^photo 
Y photo 
'^photo 


0.7265 
0.2556 
0.0179 


0.7290 
0.2514 
0.01964 


0.7050 
0.2754 
0.01948 


0.7226 
0.2584 
0.0190 


0.73507 
0.24695 
0.01798 


0.744 
0.238 
0.0182 


Rb/Rs 
T^IO*' K] 
Cb[W'' cm s"^] 


0.721 
2.142 
2.21* 


0.7115 


0.7167 

2.151 

2.19* 


0.716 
2.175 
2.22* 


0.712 
2.204 
2.25* 


0.716 

2.17 

2.22 


Tc[W' K] 
Po[100 gr cm~3] 


1.573 
1.633 


1.581 
1.559 


1.567 
1.550 


1.567 
1.545 


1.5843 
1.562 


1.569 
1.518 



Table 4: Physical inputs of the SSMs in table 3: C&:S70=Cox & Stew- 
ard (1970), MHD= Mihalas et al. (1988), CEFF=EOS of Eggleton et al. 
(1973) with the Coulombian correction added (see Christensen -Dalsgaard 
k Dappen 1992), EFF=Eggleton et al. (1973), BP92=Bahcall k Pinson- 
neault (1992), G&N93= Grevesse & Noels (1993), G91=Grevesse (1991), 
F75=Fowler et al. (1975), C&:F88=Caughlan et al. (1988). 





GGK89 


P94 


KS94 


RVCD96 


BP95 


FRANEC96 


OPACITY 


C&S70 


OPAL 


OPAL 


OPAL 


OPAL 


OPAL 


EOS 


MHD 


CEFF 


KS94 


MHD 


BP92 


OPAL 


MIXTURE 


G&N93 


G91 


G&N93 


G&N93 


G&N93 


G&N93 


CROSS SECTIONS 


F75 


BP92 


C&F88 


C&F88 


table 6 


table 6 



Table 5: Comparison among the neutrino fluxes and signals of the SSMs in 
table 3. 



P94 KS94 RVCD96 BP95 FRANEC96 



$pp [lO'^ cm 
*pep[10^ cm" 
<I>Bp [10" cm- 
*N [10^ cm" 
*0 [10^ cm" 


-■2 






59.1 


59.9 


59.4 


59.1 


59.92 


-2 


s" 




0.139 


0.138 


0.138 


0.140 


0.14 


-2 


s~ 




5.18 


4.91 


4.8 


5.15 


4.49 


-2 






0.64 


0.0611 


0.559 


0.618 


0.53 


-2 


s~ 




0.557 


0.0178 


0.481 


0.545 


0.45 


$R [10*' cm" 


■2 


s~ 




6.48 


5.83 


6.33 


6.62 


5.16 



Sca[SNU] 136.9 124.3 132.77 137.0 128 
Sci[SNU] 9.02 7.9 8^9 9.3 7.4 



Table 6: Astrophysical S'-factors [MeV barn] and their derivatives with re- 
spect to energies S" [barn] for our "best" solar model, for the "best SSM 
with diffusion" of Bahcall & Pinsonneault 1995 (BP95) and for our models. 

BP95 models (a-d) model (e) 



(our best) 




3 89x10"^^ 


3 89 X 10-25 


3 89 X 10-25 


iJ [yJ 1 Tin 


4.52x10^24 


4.52 X 10^24 


4.52 X 10-24 


5(0)33 


4.99 


5.00 


5.1 


5'(0)33 


-0.9 


-0.9 


-3.0 


5(0)34 


5.24 xlO"^ 


5.33 xlO-^ 


5.1 xlO-'' 


5'(0)34 


-3.1 xlO""* 


-3.10 xlO"'' 


-4.23 xlO-'' 


5(0) 17 


2.24 xlO-s 


2.24 xlO-s 


2.24 xlO-s 


5'(0)i7 


-3.00 xlO-5 


-3.00 xlO-5 


-3.00 xlO-5 


5(0)i2c+p 


1.45 xlO"'' 


1.40 xlO-3 


1.40 xlO-3 


5'(0)i2c+p 


2.45 xlO"-^ 


4.24 xlO-3 


4.24 xlO-3 


5(0)i3c+j, 


5.50 xlO-3 


5.50 xlO-3 


5.50 xlO-3 


5'(0)i3c+p 


1.34 xlO-2 


1.34x10-2 


1.34x10-2 


5(0)i4j\-+p 


3.29 xlO-3 


3.32x10-3 


3.32x10-3 


5'(0)i4^-+p 


-5.91 xlO--"* 


-5.91x10-3 


-5.91x10-3 


5(0)i6^r(p^^)i6o 


6.40x10-2 


6.40x10-2 


6.40x10-2 


5'(0)i5Ar(j,^^)i6o 


3.00x10-2 


3.00 xlO-2 


3.00 xlO-2 
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